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Abstract 

We analyse the asymptotic behaviour of a nonlinear mathematical model of cellular proliferation 
which describes the production of blood cells in the bone marrow. This model takes the form of a 
system of two maturity structured partial differential equations, with a retardation of the maturation 
variable and a time delay depending on this maturity. We show that the stability of this system 
depends strongly on the behaviour of the immature cells population. We obtain conditions for the 
global stability and the instability of the trivial solution. 

Keywords: Nonlinear partial differential equation, Maturity structured model. Blood production system, 
Delay depending on the maturity. Global stability, Instability. 

1 Introduction and motivation 

This paper is devoted to the analysis of a maturity structured model which involves descriptions of process 
of blood production in the bone marrow (hematopoiesis). Cell biologists recognize two main stages in the 
process of hematopoietic cells: a resting stage and a proliferating stage (see Burns and Tannock [8]). 

The resting phase, or Go-phase, is a quiescent stage in the cellular development. Resting cells mature 
but they can not divide. They can enter the proliferating phase, provided that they do not die. The 
proliferating phase is the active part of the cellular development. As soon as cells enter the proliferating 
phase, they are committed to divide, during mitosis. After division, each cell gives birth to two daughter 
cells which enter immediatly the resting phase, and complete the cycle. Proliferating cells can also die 
without ending the cycle. 

The model considered in this paper has been previously studied by Mackey and Rudnicki in 1994 pU] 
and in 1999 j21| . in the particular case when the proliferating phase duration is constant. That is, when 
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it is supposed that all cells divide exactly at the same age. Numerically, Mackey and Rey [T51[Tn], in 
1995, and Crabb et al. [51[Tn], in 1996, obtained similar results as in [20] . The model in [2D] has also been 
studied by Dyson et al [TT] in 1996 and Adimy and Pujo-Menjouet [2i Bl in 2001 and 2003, but only in 
the above-mentioned case. These authors showed that the uniqueness of the entire population depends, 
for a finite time, only on the population of small maturity cells. 

However, it is believed that, in the most general situation in hematopoiesis, all cells do not divide at 
the same age (see Bradford et al. [7]). For example, pluripotent stem cells (the less mature cells) divide 
faster than committed stem cells (the more mature cells). 

Mackey and Rey [TT], in 1993, considered a model in which the time required for a cell to divide is 
not identical between cells, and, in fact, is distributed according to a density. However, the authors made 
only a numerical analysis of their model. Dyson et al. [HITS], in 2000, also considered an equation in 
which all cells do not divide at the same age. But they considered only one phase (the proliferating one) 
which does not take into account the intermediary flux between the two phases. Adimy and Crauste [T], 
in 2003, studied a model in which the proliferating phase duration is distributed according to a density 
with compact support. The authors proved local and global stability results. 

In 2] , Adimy and Crauste developed a mathematical model of hematopoietic cells population in which 
the time spent by each cell in the proliferating phase, before mitosis, depends on its maturity at the point 
of commitment. More exactly, a cell entering the proliferating phase with a maturity m is supposed to 
divide a time r = T(m) later. This hypothesis can be found, for example, in Mitchison [55] (1971) and 
John [IS] (1981), and, to our knowledge, it has never been used, except by Adimy and Pujo-Menjouet in 
[5], where the authors considered only a linear case. The model obtained in [2] is a system of nonlinear 
first order partial differential equations, with a time delay depending on the maturity and a retardation 
of the maturation variable. The basic theory of existence, uniqueness, positivity and local stability of this 
model was investigated. 

Many cell biologists assert that the behaviour of immature cells population is an important consider- 
ation in the description of the behaviour of full cells population. The purpose of the present work is to 
analyse mathematically this phenomenon in our model. We show that, under the assumption that cells, 
in the proliferating phase, have enough time to divide, that is, T(m) is large enough, then the uniqueness 
of the entire population depends strongly, for a finite time, on the population with small maturity. This 
result allows us, for example, to describe the destruction of the cells population when the population of 
small maturity cells is affected (see Corollarv l3.ip . 

In 21J, Mackey and Rudnicki provided a criterion for global stability of their model. However, these 
authors considered only the case when the mortality rates and the rate of returning in the proliferating 
cycle are independent of the maturity variable. Thus, their criterion can not be applied directly to our 
situation. 

This paper extends some local analysis of Adimy and Crauste [2] to global results. It proves the 
connection between the global behaviour of our model and the behaviour of immature cells (to = 0). 

The paper is organised as follows. In the next section, we present the equations of our model and 
we give an integrated formulation of the problem, by using the semigroup theory. In section [31 we show 
an uniqueness result which stresses the dependence of the entire population with small maturity cells 
population. In Section (H we focus on the behaviour of the immature cells population, which satisfies 
a system of delay differential equations. We study the stability of this system by using a Lyapunov 
functionnal. In Section [5l we prove that the global stability of our model depends on its local stability 
and on the stability of the immature cells population. Finally, in Section [S] we give an instability result. 

2 Equations of the model and integrated formulation 

Let N(t, to) and P(t, m) denote, respectively, the population densities of resting and proliferating cells, 
at time t and with a maturity level to. 
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The maturity is a continuous variable which represents what composes a cell, such as proteins or other 
elements one can measure experimentally. It is supposed to range, in the two phases, from m = to m = 1. 
Cells with maturity m — are the most primitive stem cells, also called immature cells, whereas cells with 
maturity m — 1 are ready to enter the bloodstream, they have reached the end of their development. 

In the two phases, cells mature with a velocity V{m), which is assumed to be continuously differentiable 
on [0, 1], positive on (0, 1] and such that V{0) = and 



/' 

Jo 



ds 

— — = +00, for me (0,11. (1) 
V{s) 



Since j™^ vfs)^ ^ith mi < m2, is the time required for a cell with maturity mi to reach the maturity m2, 
then Condition ([T]) means that a cell with very small maturity needs a long time to become mature. 
For example. Condition ^ is satisfied if 

V{m) ~ am^, with a > and p >1. 

m— »0 

In the resting phase, cells can die at a rate 5 = 5{m) and can also be introduced in the proliferating 
phase with a rate (3. In the proliferating phase, cells can also die, by apoptosis (a programmed cell death), 
at a rate 7 = 7(m). The functions 5 and 7 are supposed to be continuous and nonnegative on [0, 1]. 
The rate (3 of re-entry in the proliferating phase is supposed to depend on cells maturity and on the 
resting population density (see Sachs [13]), that is, /3 = /3(m, N{t, m)). The mapping (3 is supposed to be 
continuous and positive. 

Proliferating cells are committed to undergo mitosis a time r after their entrance in this phase. We 
assume that r depends on the maturity of the cell when it enters the proliferating phase, that means, if a 
cell enters the proliferating phase with a maturity m, then it will divide a time r — T[m) later. 
The function r is supposed to be positive, continuous on [0,1], continuously differentiable on (0,1] and 
such that 

^'M + 77^ > 0' for e (0> !]■ (2) 
V[m) 

One can notice that this condition is always satisfied in a neighborhood of the origin, because V{Q) — 0, 
and is satisfied if we assume, for example, that r is increasing (which describes the fact that the less 
mature cells divide faster than more mature cells). 
Under Condition ([2]), if m G (0, 1] is given, then the mapping 

is continuous and strictly decreasing from (0,m] into [—r(m), +00). Hence, we can define a function 
8: (0,1] ^(0,1], by 

ds 

/ T77T = ^(0(H), for me (0,1]. 

The quantity 0(m) represents the maturity of a cell at the point of commitment when this cell divides 
at a maturity level m. The function G is continuously differentiable and strictly increasing on (0, 1] and 
satisfies 

lim 0(m) — and < 0(m) < m, for m e (0, 1]. 

m— ^0 

If we consider the characteristic curves x ■ (~oo,0] x [0, 1] [0, 1], solutions of the ordinary differential 
equation 

{dx 
— (s,m) = y(x(s,m)), s < and m e [0, 1], 
ds 
X(0,m) = m, 
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then, it is easy to check that, for m e [0, 1], ©(m) is the unique solution of the equation 

X = x{-Tix),m). (3) 

The characteristic curves x(S)TO) represent the evolution of the cell maturity to reach a maturity m at 
time from a time s < 0. They satisfy x{s, 0) = and x(s, m) S (0, 1] for s < and m G (0, 1]. Moreover, 
we can verify that the characteristic curves are given by 

X{s,m) = h~^{h{m)e^), for s < and m G [0, 1], (4) 

where the continuous function h : [0, 1] [0, 1] is defined by 

Mm) = | -P(-Iv^)' forme (0,1], 
[ 0, form = 0. 

Since /), is increasing, the two fmictions s t—^ x(s,to) and m ^ x(,s. m) arc also increasing. 

At the end of the proliferating phase, a cell with a maturity m divides into two daughter cells with 
maturity g{m). We assume that g : [0, 1] [0, 1] is a continuous and strictly increasing function, continu- 
ously diffcrcntiablc; on [0, 1) and such that g{m) < m for m G [0, 1]. We also assume, for technical reason 
and without loss of generality, that 

lim g'{m) = +oo. 

m— »1 

Then we can set 

fif"^(m) = 1, for m> 5(1). 

This means that the function g~^ : [0, 1] [0, 1] is continuously differentiable and satisfies {g~^y{m) = 0, 
for m > g{l). We set 

A{m) = e{g-\m)), form €[0,1]. 

The quantity A(m) is the maturity of a mother cell at the point of commitment, when the daughter cells 
have a maturity rn at birth. The function A : [0, 1] — > [0, 1] is continuous and continuously differentiable 
on (0, 1]. It satisfies A(0) = 0, A is strictly increasing on (0,^(1)), with 6(to) < A(m), and A(m) = 6(1) 
for mG [5(1), 1]. 

At time t = 0, the resting and proliferating populations are given by 

N{0,m) =Jl{m), (5) 

and 

l-r{e{m)) 

P(0,m) = r(m) := / T{m,a)da, (6) 

./o 

where T[m,a) is the density of cells with maturity m, at time t = 0, which have spent a time a in the 
proliferating phase, or, equivalently, with age o. The functions Jb and T are supposed to be continuous on 
their domains. 

We define the sets 

f2:= [0,1] X [0,w], 

where Tmax ■= max„g[o,i] t(to) > 0, 

:= |(rn,t) G ; 0<t< r(A(m))|, 

and 

:= |(TO,t) en ; 0<t< r(e(m))|. 
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(7) 



Then, the population densities N{t,m) and P{t,m) satisfy, for rn E [0,1] and t > 0, the foUowing 
equations, 

^N(t,m) + -^{V{m)N{t,m)) = -(^S{m) + f3{m, N{t,m))'jN{t,m) 
' 2e(t,TO)r(x(-t,ff-^(TO)),r(A(m))-i), if(m,i)el7A, 
2^{T{A{m)),m)p(A{m), N{t - T(A(m)), A(m))) x 



and 



(8) 



N{t^T{A{m)),A{m)), if(m,i)^riA, 
d d 

-^P{t, m) + —{V{m)P{t, to)) = ~j{m)P{t, m) + /3(to, N{t, m))N{t, to) 

L/C CJift 

' TT{t,m)T(^x{ -t,m),T{Qim))-ty if (to,*) G fie, 

< 7r(r(e(m)),m)/3(e(m),iV(t-T(e( to)), 0(m)) 



iV(t-T(e(TO)),e(m)), if(TO,<)^17e, 
where the mappings ^ ; J^a ^ [0, +oo) and tt : f2e ^ [0, +oo) are continuous and satisfy 

to) = if TO > g(l), 

because, from the definition of g, a daughter cell can not have a maturity greater than 17(1). 

In Equation (O, the first term in the right hand side accounts for cellular loss, through cells death (6) 
and introduction in the proliferating phase (/3) . The second term describes the contribution of proliferating 
cells, one generation time ago. In a first time, cells can only proceed from cells initially in the proliferating 
phase (r). Then, after one generation time, all cells have divided and the contribution can only comes 
from resting cells which have been introduced in the proliferating phase one generation time ago. 
The factor 2 always accounts for mitosis. The quantity ^{t,m) is for the rate of surviving cells. 
In Equation ([S]), the first term in the right hand side also accounts for cellular loss, whereas the second 
term is for the contribution of the resting phase. The third term describes the same situation as in 
Equation ([7]), however, in this case, cells leave the proliferating phase to the resting one. The quantity 
7r(i, to) is also for the rate of surviving cells. 

We can observe two different behaviours of the rates of surviving cells, in the two phases. In a first time, 
they depend on time and maturity, and after a certain time, they only depend on the maturity variable. 
When the process of production of blood cells has just begun, the only cells which divide come from 
the initial proliferating phase population. But after one cellular cycle, that means when t > t(A(to)) 
(respectively, t > t(0(to))), the amount of cells only comes from resting cells (respectively, proliferating 
cells) which have been introduced in the proliferating phase (respectively, resting phase) one generation 
time ago. Consequently, we take into account the duration of the cell cycle, and not the present time. 
Equations ((T]) and ([5|) are derived, after integration, from an age-maturity structured model, presented by 
the authors in [2]. In fact, the rates ^ and tt are explicitly given (see [I]) by 



at,m) = (.g-i)'(TO)exp| - / [j{x{-s, 9'\m))) + V {xis, g-\m)))^ 



ds 
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and ^ 

n{t, to) = a{m) exp J (i{x{~^t''t^)) + ^' {x{~s,m))^ds 

with a : [0, 1] [0, +oo) a positive and continuous function, such that a(0) = 1. 
In the foUowing, to simphfy the notations, we will denote by ^ and 7f the quantities 

l{m) = ^(r(A(TO)),m), 

and 

Tr{m) — 7r(r(0(TO)), to) . 

One can remark that the solutions of Equation ([7]) do not depend on the solutions of Equation ^ , whereas 
the converse is not true. 

Before we study the asymptotic behaviour of the solutions of Problem (O-®, we establish an integrated 
formulation of this problem. We first extend N by setting 

N(t, to) = Jl{m), for t e [~T,nax,0] and to G [0, 1]. (9) 

One can remark that this extension does not influence the system. 

We also define two mappings, F : [0, +oo) x [0, 1] x M ^ R and G : [0, +oo) x [0, 1] x M ^ R, by 

[ 2at,m)r(x{-t,g-\m)),T{A{m))-t), if(m,i)er!A, 
F{t,m,x)^< (10) 

[ 2l{m)l3{A{m),x)x, if (m, t) i ^a, 

and 

{7r{t,m)T(x{-t,m),T{e(m)) -t), if {m,t) G f^e, 
(11) 
Tf{m)(3(Q{m),x)x, if {m,t) ^ Hq. 

We denote by C[0, 1] the space of continuous functions on [0, 1], endowed with the supremum norm ||.||, 
defined by 

||t;|| = sup |w(to)|, for t; G C[0, 1], 

me [0,1] 

and we consider the unbounded closed linear operator A : D{A) C C[0, 1] — » C[0, 1] defined by 
D{A) = iue C[0, 1] ; u differentiable on (0, 1], u' G C(0, 1], lim V{x)u{x) ^ o] 



and 

-{S{x) + V'{x))u{x) - V{x)u'{x), if a; G (0, 1], 
-((5(0) + y'(0))u(0), ifa; = 0. 



Au{x) 



Proposition 2.1. The operator A is the infinitesimal generator of the strongly continuous semigroup 
{T{t))t>o defined on C[0, 1] by 

{T{t)4>){x) ^ K{t,x)^p{xi-t,x)), for4>eC[0,l], t > and x e [0,1], 

where ^ 

K{t, x) = exp |- ^ {S{xi~s, x)) + V'ixi-s, a;)))ds 

Proof. The proof is similar to the proof of Proposition 2.4 in Dyson et al. [TT]. □ 
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Now, by using the variation of constants formula associated to the Co-semigroup (T(t))t>o, we can write 
an integrated formulation of Problem (O-®. 

Let C{Q) be the space of continuous functions on Q, endowed with the norm 

||T||n:= sup |T(m,a)|, for T e C(r!). 

Let JI G C[0, 1] and F G C(fl). An integrated solution of Problem (H])-® is a continuous solution of the 
system 

N{t,m) = K{t,m)Ji{x{-t,m)) 



K{t - s,m)/3(x(-(t - s),m),N{s,x{-{t - s),m))^N {s,x{-{t - s),m))ds 
+ j^K(t- s, m)F(s, x{-{t - s), m), N{s ~ r(A(x(-(t - s), m))), A(x(-(t - s), m)))) 

and 

P(i,m) = H{t,m)T{x{~t,m)) 

+ I H{t-s,m)l3(x{-{t-s),m),N{s,x{~{t~s),m))^N{s,x{-{t-s),m))ds 



ds. 



H{t - s, m)G(s, x(-(t - s),m),N{s - T(e(x(-(t - s),m))), e{x{-(t - s), m))))ds, 
for t >0 and to G [0, 1], where F and G are given by lfTO|) and ([TTjl . F is given by ([6]) and 



H{t,m) := exp (j{x{-s,m)) + V'{x{-s,m))y 



ds> , for i > and to G [0, 1]. 



(12) 



(13) 



We can easily prove (see [2]), under the assumptions that the function x I3{m, x) is uniformly bounded 
and the function x x/3{m,x) is locally Lipschitz continuous for all m G [0, 1], that Problem (fT^ - (fT5|) 
has a unique continuous global solution (iV'^''", P^'^), for initial conditions (/I, F) G C[0, 1] x C{Q). 



3 A uniqueness result 

In this section, we establish more than uniqueness. Indeed, we show a result which stresses, for a finite 
time, the dependence of the entire population with the small maturity cells population. It has been shown 
for the first time by Dyson et al. [11 , for a model with a constant delay. We will see that this result is 
important in order to obtain the asymptotic behaviour of the solutions of (fT2|l - (ll3p . 
We first assume that 

A(m) < TO, for all m G (0, 1]. (14) 

This condition is equivalent to 

ds 

t(A(to)) > / -— , for TOG (0,1]. (15) 

Jrn y[-V 

This equivalence is immediate when one notices that, from ([3]), 

A{m)=x{-T{A{m)),g~\m)) = h-\h{g-\m))e-''^'^^^"''>^) . 
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Since the quantity represents the time required for a cell with maturity m, at birth, to reach 

the maturity of its mother at the cytokinesis point (the point of division), Condition (jlSp means that, in 
the proliferating phase, cells have enough time to reach the maturity of their mother. 
Condition (|14p implies in particular that 

6(1) :=A(5(1))<5(1). 

From now on, and throughout this section, we assume that the function x i— > f3{m, x) is uniformly bounded, 
the function x i— > x(3{m, x) is locally Lipschitz continuous for all m G [0, 1], and that Condition (HH) holds. 
For b e (0, 1] and -tjj G C[0, 1], we define ||.||b as follows 

llV^llb := sup |V'(™)|- 

me[Q,b] 

We first show the following proposition. 

Proposition 3.1. Let ]Ii,]l2 G C[0, 1] and ri,r2 € C(ri). // there exists < 6 < 1 such that 

'Pi{m) = 'p2{ra) and Ti(m,a) = T2(m,a), (16) 

for m G [0, b] and a £ [0, Tmax], then, 

m^''^^{t,m) = Nf^^-^^{t,m), fort>Oandme[0,g{b)]. (17) 

Proof. We suppose that there exists b € (0, 1) such that (fTB]) holds. Let T > be given, and let t G (0, T] 
and m G [0,5(6)] be fixed. Since h is increasing, it follows from ^ that 

x(— m) < m < g{b) < b. 

Then 

Let s G [0,t]. Since g^^ is increasing, then 

x{-s,g-\x{-{t-s),m))) < g'^ {x{-{t - s),m)) <g-\m)<b. 
Moreover, if < s < t(A(x(— (i — s), m))) , then 

T(A(x(-(t-s),m))) - s G [0,r™,J. 

Thus, we have 

ri(x( - s,g-^{x{-{t - s),m))),r(A(x(-(t - s),m))) - s) 

= r2 (x( - g-^ {xi-it - s),m))), T(A(x(-(t - s),m))) - s) . 
Since the solutions N^^-'^'^{t,m) and N^^'^^{t,m) of Equation (fT2|) are continuous and satisfy 

m^^^^{0,m) = 7V^2>r2(o^m)^ for m G [0,6], 
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then, by using the locally Lipshitz continuous property of the function x i-^ xP(m,x), we can write 
\N'f^^'^^{t,m) - ^^•^^{t,m)\ 

<KL I \m^'^'{s,xi-it - s),m)) - m^'^''{s,x{-{t - s),m))\ds 
Jo 

+2KLm f |iV^^-r, _ r(A(x(-(t - s), m))), A(x(-(i - s), m))) 

-^2X2^, _ r(A(x(-(i - s),m))), A{x{-{t - s), m)))|ds, 

JO 

+2KLm\ f llTV^-i^n^^ - T(A(x(-(i - 5),m))), .) - 7V^^-i^^(s - T(A(x(-(t - s),m))), 
Jo 

for T > small enough, where L is a Lipschitz constant of the function x i— s- x(3{m^ x) and K is defined 

by 

K{s, m) < K, for s £ [0, T] and m £ [0, 1]. 

The extension given by ^ allows to give sense to the integral terms in the above inequality. 

Let e e [-TynaxA be given. If t + 61 < 0, then m^'^^[t + O^m) = m^'^-^{t + O^m). If t + 6I > 0, then 

\mi^^^{t + e,m) - m-^^^^{t + e,m)\ 

Jo 

+2KLm\ / llTV^-i^n^ - r(A(x(-(t + 6 - s),m))), .) - N^-^^^is - T(A(x(-(t + - s),m))), 
Jo 

<i^L(l + 2lje||) r sup \\mi'^^is + 9,.)-m-^^-is + 9,.)\\g(^ds. 

JO 0e[-T„„,,o] 

It follows that 

sup \\N^^^^^ {t + 0, .) - m-^^-{t + 0, 

e6[-r„„,,0] ^ 

<ifL(l + 2||e||) / sup ||iVA^i^r,(^^^^ )_^7i„r.(^_^^^_)|j^^^^^^^ 

Jo ee[-r,„a^,o] 

By using the Gronwall's Inequality, we obtain 

sup \\m^^^^{t + e,.)-NT'-^''-{t + 6,.)\\g^,)=0. 

6'e[-r„„^,0] 

In particular, 

\\NT'^''''{t,.)-N^-^^-{t,.)\\g(b)^Q. forte (0,T]. 
By steps, this result holds for all T > 0, therefore (flTl) is satisfied and the proof is complete. □ 



9 



M. Adimy, F. Crauste and L. Pujo-Menjouet 



Stability of a cellular proliferation model 



Now, let < 6 < g{l) be fixed and consider the sequence (6„)„gN defined by 

A-i(6n), if 6„ e [0,6(1)), 
5(1), if6„e [e(l),5(l) 



bo = b and 6„+i 



(18) 



The sequence (6„)„gN represents the transmission of the maturity between two successive generations, n 
and n + 1. The following result is immediate. 

Lemma 3.1. If Q <b < 6(1) := lS,{g{l)), then there exists iV e N such that b^ < 6(1) < bN+i < .9(1)- 

We give now a first result, which emphasizes the strong link between the process of production of cells 
and the population of stem cells. A similar result has been proved by Adimy and Pujo-Menjouet 5 in 
the linear case. 

Theorem 3.1. Let /ii:M2 ^ ^"[0,1] and ri,r2 G C{fl). If there exists < 6 < 1 such that II 16)) holds, 
then, there exists t > such that 

m^'^'it,m) = m^^^^t,m), 
for m e [0,5(1)] and t >t, where t can be chosen to be 

"M5(l)) 



t = In 



+ {N + 2)r„ 



H9ib))_ 

and N N is given by Lemma \8.1\ for b = g{b). Furthermore, 



(19) 



iV^i^^i(i,m) ^ Ar^2^^^(i,m), 
for m G [.g(l), 1] andt>t + T„,ax - In (/^(.9(1))) = {N + i)T„,ax - In {h{g{b))) . 
Proof. Let b = g{b). Since g is increasing, then b < 5(1). Proposition 13. II implies that 

A^'^i '^1(^,771) = N'f^^^^''{t,m), for t > and m e [0,6]. 
Let us reconsider the sequence (6n)nGN, given by p^ . and let us consider the sequence (i„)„gN defined by 

h{bn+i) 




h{bn 



(20) 



Then, 



tn = In 



Hbr^) 

Hg{b)) 



UTr, 



The sequence (fen)neN is increasing. Then, the sequence (in)neN is also increasing. We are going to prove, 
by induction, the following result 

(i?„) : iV'^i'^i(t,m) = A^^2J2(^^^)^ for i > i„ and m e [0,5„]. 



First, [Hq) is true, from Proposition 13. II 

Let suppose that (-ff„) is true for n e N. Let t > t„+i and m e [0, 6,1+1]. Then, from (PO)) . 

tn+l ^ tji -\- T^iQ^x ^ '^raax- 
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ds 



+2 



Since Equation p2p is autonomous, its solutions can be reformulated, for t > as follows 

iV^^'^*(i,m) = K{t - tn - T^ax^m)^'^'^' {tn+Tmax,x{ ^ (t - tn - Tmax) , m)^ 

K{t-s,m)p(^x{-it- s),m) , N^^'^^ (5, x( - (t - s), m)) ^ iV^^-r, (^^^ ^( _ _ 

K{t - s, m)e(x( -it- s), m))iV^^-r, _ T(A(x(-(t - s),m))),A{x{-{t - s), m))) x 

(3 (^A{x{-{t - s), m)), iV^^-r, _ r(A(x(-(t - s), m))) , A(x(-(t - s), m))) ^ ds, 
for i = 1,2. Remark that, from (jl]), 

and from (PU)) . 

~ h{bn+i)' 

Then, we deduce that 

xi - it-tn-Tmax),m) < h^^ ( h{m) ,^}^"W 

Hence, (Hn) implies 

Nf^^'^^ (^tn + Tmax,x{ - (t - t^ - Tmax) , m)) = N'^2,r2 (^^ Tmax,x{ - [t - tn - T^ax) , w) j . 

Furthermore, for tn + Tmax < s < t, we have 

s - r(A(x(-(t - s), m))) > (t„ + r„,„^) - r(A(x(-(i - s), m)) > i„, 

and 

A(x(-(^ - s),m)) < A(m) < A(6„+i) = &„. 

Consequently, 

7VM„ri (,_^(A(x(-(i-s), m)), A(x(-(i-s), m))) = 7VM..r. (s-r(A(x(-(i-s), m)), A(x(-(t-s), m)))) . 
Then, we obtain that 



t 

K{t — s, m) 



P{x{~{t - s), m), iV^-^i _ s), m)))iV^-ri _ 

■P{xi-it - s),m),m-'^- (s, x(-(i - s),m))^N^^^^- {s, x(-(t - -s), m)) 
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and, by using the Gronwall's incquahty, we deduce that {H„+i) is true. Consequently, (Hn) is true for 
n G N. 

In particular, {Hn) holds for n — N + 2, where N is given by Lemma [3.11 with b = g{b). Since 6Ar+i e 
[9(1), 5(1)], then bN+2 5(1)- We deduce that 

7V'^i'^i(i,m) = N'^^^^^{t,m) for m G [0,g(l)] and t > t, (21) 

where t — t is given by UH) . 

Finally, take m G [5(1), 1] and t >t + T,nax- We can write, for i — 1,2, 

N't^'-'^^it, m)^K{t-t- T„,ax,m)m^''^^{t + r„„^, x(-(t -t- T,nax),m)) 

Kit ~ s, m)[i(^x{ -{t- s),m),m^'^^ {s,x{ - (t ~ s), m))^iV^^-r, (^^^ ^(-(t- s), m))ds 
K{t - s, m)e(x( - (t - s), m))iV^.-r- (s - T(A(x(-(t - s) , m))) , A{x{- (t - ,s), m))) x 

p(^A{x{-{t - s), m)), TV^^-r. - T(A(x(-(t - s), m))) , A(x(-(t - s), m)))^ 

Let < + Tmax < s < Then, 

s-r(A(x(-(t-s),m))) > (t + w) - r(A(x(-(i - s), m)) >I 



-2 



t 

t + T„ 



ds. 



Consequently, if 
then, 



x{-{t-s),m) <g{l), 
A(x(-(i-s),m))<A(5(l))<.9(l), 



and (PT|) implies that 

iVMiTi (s-r(A(x(-(t-s), m))) , A{x{-{t-s), m))) = iV^^-^, (^s-T{A{x{-{t-s), m))) , A(x(-(i-s), m))) . 
On the other hand, if 

x{ -it-s),m) >5(1), 

then, from the definition of ^, we have 

C(x(-(i-s),m)) -0. 

Furthermore, by remarking that hi(h(m)) < 0, for all m G (0, 1], then we deduce, for m G [5(1), 1] and 
t>t + Tmax - ln(ft.(g(l))), that 

X(- (i-t- w),m) = /i-i(/i(TO)e-(*-*-"— 

< /i-i(/i(m)/i(5(l))), 

< /z-i(%(l)))=g(l). 



Hence, 



^''1^^' + Tmax.x[ -it-t- Tmax). m) ) = A^^^^^^ (t + r,„ax, x( " (^ " * " T^ax). m) ) . 
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Using once again the Gronwall's inequality, we conclude that 

N^^'^^{t,m)^NT'-^^^it,m), for m S [r7(l),l] and t > I + w - ln(/i(5(l))). 
This completes the proof. □ 
Corollary 3.1. Let Jl e C[0, 1] and T e C{n). If there exists <b <1 such that 

Jl{m) = and T{m, a) — 0, for m G [0, b] and a £ [0, Tmax]^ 

then, 

N'f^'^{t,m) = 0, for m G [0, 1] and t > {N + 3)Tmax - hi {h{g{b))), 
where N ieN is given by Lemma \3.1\ for b — g{b). 

This result stresses the dependence of the production of cells with the population of stem cells. In 
particular, if the stem cells population is defective in the initial stage, then the entire population is doomed 
to extinction in a finite time. This situation describes what usually happens with the aplastic anemia, a 
disease which yields to injury or destruction of pluripotential stem cells. 

In the next corollary, we show that the proliferating population depends also strongly on the stem cells 
population. 

Corollary 3.2. Let 'Jli,'p2 C'fO, 1] and Ti,T2 E C{Q). If there exists < b < 1 such that jT^j holds, 
then 

P^i'^i(t,m) = F^2'^^(i,m), forme [0,1] and t > {N + 3)Tmax - '^n {h{g{b))) , 
where N £ N is given by Lemma \3.1[ for b — g{b). 

Proof. The proof is immediate by using Theorem 13. 11 Equation (113p and a method of steps. □ 

4 Behaviour of the immature cells population 

In this section, we investigate the behaviour of the immature cells population, that means, the population 
of cells with maturity to = 0. 

Let Jl € C[0, 1] and F G C{il) be fixed. Let us consider the continuous solutions Af^''"(t,TO) and 
pT^-^{t,m) of Problem HI])-!!!]). We set x{t) = ^■^{t.'d) and y{t) = P'^^■^{t,Q), for all t > 0. Then, 
{x{t),y{t)) is solution of the system 



xit) = e-P'-piO)- / e-P'-'-''>(3iO,x{s))x{s)ds 
t ° 

2 I e"''(*-'')^(s,0)r(0,r-s)ds, foriG[0,r] 



2 / e-''(*-'')C(s,0)r(0,r-s)ds + 2^(0) / e-''('"")/3(0, x(s - r))a;(s - r)rfs, for r < i, 

(22) 



and 

y{t) = e-''*r(0)+ / e-'"^^-'^l3{0,x{s))x{s)ds 







g-^(t-.)^(g^O)r(0,r-s)ds, for<G [0,r], 

''(*-'')7r(s, 0)r(0, r - s)ds + 7f(0) / e"''(*~'*)/3(0, x{s - r))x{s - r)ds, for r < i. 



e 







(23) 
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where p := S{0) + V'{0), rj 7(0) + V'{Q) and r := t(0) > 0. 
Let us recall that ^(0) = kir, 0), 7f(0) = 7r(r, 0) and 

r(0) = I r(0, a)da. 
Jo 

Then, we easily deduce that System ([22l) -(|23 |) is equivalent to the system 



dx^^^ _ ^ ata , / 2C(t,0)r(0,r - t), for < t < r, 

, ioir <t, (24) 



-it) . -(p + /^(o,^W)).(t) + |;|--^^^^^^^ 

x(0) = 71(0), 



^(t) = -vy(t) + mx(t)Ut)-h^''^^^^^''-'^' forO<i<r, 
dt^' ^Vy[t) + P[^,H^))^[V \ 7r{0)p{0,x{t~r))x{t-r), for r < (25) 

y(0) = r(o). 

Of course, at i = r, the derivatives in (|24p and (P5|) represent the right-hand side and the left-hand side 
derivatives. 

First, consider the system, for t e [0,r], 

^(t) = -(p-H/3(O,0(O))0(i) + 2e(t,O)r(O,r-i), 
-^{t) = -Tim+P{0,<P{tmt)-7^{t,0)r{0,r-t), 
with 

0(0) = ^(0), 

V'(o) = r(o). ^ ^ 

It is obvious that, under the assumptions that the function x 1-^ /3(0, x) is bounded and the function 
X 1-^ x(3{0,x) is locally Lipschitz continuous. Problem (P51l - (|27p has a unique solution (^(f>(t),^(t)) , for 
t e [0, r]. Remark that ip{t) is explicitly given by 

V;(t)=e-''*/ r{0, a)da+ e-''(*-'*)/3(0, 0(s))<?!)(s)ds, forie[0,r]. (28) 
^0 Jo 

Moreover, if /l(0) > and r(0, •) > 0, then (p{t) and ip{t) are nonnegative. 
Hence, for t > r, Problem reduces to the delay differential system 

^(t) = -{p + /?(0, x{t)))x{t) + 2?(0)/3(0, xit ~ r))x{t ~ r), (29) 
^(t) = -wit) + P{0,xit))xit) -W{0)PiO,xit - r))x{t - r), (30) 

with, for t e [0,r], 

x{t) = 0(t), 



y{t) V(i). 
As ip{t), for t G [0, r], ?/(t) is explicitly given, for t > r, by 



(31) 



y{t)= / e-''(*-")/3(0,x(s))a;(s)ds. (32) 
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Proposition 4.1. Assume that the function x i~> (3{0,x) is bounded and the function x >— > xP{0,x) is 
locally Lipshitz continuous. Lefp G C[0, 1] and T G C{Q) be given. Then, Problem i29\) - ^^) has a unique 
solution {x'^ {t) ^ y^ [t)) , defined for t > 0, where (^(j){t),ip(t)^ is the unique solution of h26\) - \27^ . Moreover, 
{x'^{t),y'^ (t)) has a continuous derivative at t = r if and only if 

r(0,0)=/3(0,/l(0))7l(0). (33) 
Furthermore, if]l{Q) > and r(0, •) > 0, then x'^{t) and y'^{t) are nonnegative. 

Proof. Existence, uniqueness and regularity of solutions of Problem (j29|) - (l31l) come from Hale and Verduyn 
Lunel m]. The positivity of x'^{t) is easily obtained by steps. Moreover, by using ([5^ . we deduce 
immediatly the positivity of y^'{t). □ 

In the sequel, we will consider, for a biological reason, only nonnegative solutions of Problem ((29 |) -([30 |l . 

Lemma 4.1. // linii^+oo a;'^(t) = C exists, then \mit->+ooy^{t) exists and is equal to 

i(l-e-'"-)/3(0,C)C, */ry>0, 
_ r(3{0,C)C, tfv^O- 
Proof. We assume that 

lim x'f'it) = C. 

t — * + CxD 

By using ((5^ . we obtain that 

y'^{t)= I e-''^'l3{0,x'^{t- s))x'^{t- s)ds, for t > r. 



Then, we easily conclude that 



lim^/(i)=(^^ e-'"'dsj/3(0,C)C. 



t — *-|-oo 

This ends the proof. □ 

Lemma I^TT] implies that, in order to study the stability of the solutions of Problem ((^^ - ((Sn)) . we only 
need to concentrate on the stability of the solutions of the delay differential equation (1^^ . 
In [16], Mackey has proposed that the hmction /3(0, •) is a Hill function, defined by 

nn 

where /3o and 9 are two positive constants and n > 1. This function is used to describe, from a reasonable 
biological point of view, the fact that the rate of re-entry in the proliferating compartment is a decreasing 
function of the total number of resting cells. 

We recall that the function x <—>■ P{0,x) is supposed to be continuous and positive. From now on, we also 
suppose that x i-^ P{0,x) is decreasing on [0, -|-cx)) and satisfies 

lim (3{0,x) = 0. (35) 

These assumptions have been done, for the first time, by Mackey [T^ in 1978 and have been used by 
Mackey and Rudnicki [20^ in 1994. 

Before studying the stability of Problem we recall a non-trivial property of the solutions of (^5)1 . 

The result in Proposition l4.2l has been proved for a similar equation by Mackey and Rudnicki [2^, in 1994. 



15 



M. Adimy, F. Crauste and L. Pujo-Menjouet 



Stability of a cellular proliferation model 



Proposition 4.2. Assume that p> 0. Then, every solution of Equation is hounded. 

One can notice that, if p = 0, then Equation (j29p may have unbounded solutions. A counter example 
is given in the next proposition. 

Proposition 4.3. Assume that p — and that there exists a; > such that the function x i— > x(3{0,x) is 
decreasing on \x, +oo). Let Ji G C[0, 1] and T G C(r2) he such that fiSS\] holds, 7l(0) > x and 

2C(t,0)r(0,r-i) > r(0,0), /orie[0,r]. (36) 

Then, the solution of Equation h29\) is unbounded. 

Proof. Let consider the solution x{t) of the problem 

x'{t) ^ 2^{0)(3{0,x{t-r))x{t-r) - (3{0,x{t))x{t), for t > r, 
x{t) = forO<t<r. 

First, one can notice that, if lim(_^_|_oo x{t) exists and is equal to C, then C — 0. 
By contradiction, if we suppose that C > 0, then we obtain that 

lim x'{t) = (2^(0) - l)/3(0, C)C > 0, 

t—>+oo 

because Condition p6|) implies that 2^(0) > 1. This contradicts the fact that x{t) converges. Then C — 0. 
Secondly, let Jl e C[0, 1] and T e C(f7) be such that ^ and ^ hold, and /l(0) > x. The solution (j){t) 
of the problem 

= 2at,0)r{0,r~t)~ mnmit): fori>r, 

m = 71(0), 

satisfies 

0'(o) = 2e(o,o)r(o,r)-r(o,o)>o. 

Consequently, there exists e e (0,r] such that (j)'{t) > 0, for t G [0, e). 
Hence, 0(0) < 0(e). Then, 

0'(e) > r(O,O)-/3(O,0(e))0(e), 

> /3(O,7l(O))7l(O)-/?(O,0(e))0(e), 

> 0. 

By steps, we conclude that (t)'{t) > 0, for t G [0,r]. 
By the same way, we obtain that 

x'{r) > /3(o, m)m - m > o. 

By using the same reasonning, we show that x'{t) > 0, for t > r. Hence, x is unbounded and the proof is 
complete. □ 

Remark 1. Even if the trivial solution of h29\) is unstable, the trivial solution of ViO\) may be stable. For 
example, if\m\t—,+ooX'^{t) = +oo, then, by using hSS]) and !i35\) . we easily obtain that\m\t^-^-ooy^{t) = 0. 

The assumption on the function x i-^ x(3{Q, x), in Proposition 14.31 holds, for example, if /? is given by 
with n > 1. In this case, the function x ^ xP{0, x) is decreasing for x >x = 9/{n ~ 1)^/". 
We determine, in the next theorem, the global stability area of the trivial solution of Equation 

Theorem 4.1. The trivial solution of Equation \29fl is globally stable if and only if 

(2C(0)-l)/3(0,0)<p. 
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Proof. First, we assume that (2^(0) — l)/3(0, 0) < p. We are going to show the global stability by using a 
Lyapunov functional. 

We denote by C+ the subset of C[0, r] containing nonnegative functions. We set 



f{x)=xl3{0,^) 
j:{x) = / f{s)ds, for x > 0. 



and 

We define the mapping J : ^ M by 

J(0) = T{(t>{r)) +e(0) r f{H<^))da, 
Jo 

Then 
Since 



for 6eC^ 



JW = -^{r)fi^{r))+m{fWr)) ~ f\m))- 



-f^ir) = ~p<\>[r) ~ /(0(r)) + 2C(0)/(<^(0)) , 



then 



J(0) 



p/3(0, mWir) - /2(0(r)) + m{j\4>ir)) + 2/(0(r))/(0(O)) - f\m) 
p + /3(0, 0(r))) /3(0, (r) + 2mfHHr)) ~ ?(0) f /(^W) - /(^(O)) 



J{4>) < -(p-{2m - l)/3(O,0(r)))/3(O,0(r))02(^) 



Hence, 



Since (2^(0) — l)/?(0, 0) < p and the function x /3(0,a;) is decreasing and positive on R+, then the 
function 

X{u) = (p- (2^(0) - l)/3(0,u))/3(0,u)u2 



is nonnegative on and A(u) = if and only if u = 0. Consequently, every solution of Equation 
with (p G C+ tends to zero as t tends to +oo. 
Now, if we assume that p < (2^(0) — l)/3(0, 0), then, immediatly, 

-ip + mO)) > -2C(0)/3(0,0). 

Hence, by using Bellman and Cooke ([6], Theorem 13.8), we obtain that the trivial solution of (|29|) is 
unstable. □ 

Remark that, by using Lemma [4.11 if (2(^(0) — l)/3(0, 0) < p, then the trivial solution of ([50]) is also 
globally stable. 

We are going to use the results of Theorem 14.11 in the next sections to obtain global stability and 
instability for the solutions of Problem (fT^ -(fTc 
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5 Global stability for the maturity structured model 

In this section, we establish a result of global stability for Problem P^ -(|13 p which stresses the influence 
of immature cells on the total population. First, we recall some definitions. 

Definition 5.1. The trivial solution of Problem US^} - il3]) is locally stable if, for all £ > 0, there exist 
ly > and T > such that, ifJlG C[0, 1] and T G C(f2) satisfy 

\\J1\\ < v and \\T\\n < v, 

then 

|liV^'^(t,.)|| < £ and \\P'^'■^{t,.)\\<e, fort>T. (37) 

The trivial solution of Problem hV^) - ]iT3i) is globally stable if, for all £ > 0, there exists T > such that 
^ holds. 

Throughout this section, we assume that the function x i— > /9(to, x) is uniformly bounded and that the 
function x i— > x[3{m, x) is locally Lipschitz continuous for all m G [0, 1]. In the next theorem, we show our 
main result, which makes the link between the global stability of the trivial solution of Problem (fT2|) - (fT3|) 
and the stability of the immature cells population. 

Theorem 5.1. Assume that Condition {1^-^ holds. Let us suppose that the trivial solution of Problem 
W0j) - n^) is locally stable. Then this solution is globally stable on the set 

17gs = |(M,r) G C[0,1] X C(f7) : lim Ar^'^(t,0)= lim 0) = o|. 

Proof. We first show that, if the trivial solution of Equation is locally stable, then it is globally stable 
on the set 

Ojv = ((7i,r) G C[0,1] X C(f}) : lim iV'''r(t,0) = 0|. 

Let us suppose that the trivial solution of (fT2|) is locally stable. Then, for all JI G C[0, 1], F G C{VL) and 
£ > 0, there exist v > and T > 0, such that, if 

||7I|| < V and ||F||a < v, 

then, 

|7V'^'^(i,m)| < £, fort >randmG [0,1]. (38) 
Let £ > be given and let (/I, F) G fijv- Then limt^+oo N^'^{t, 0) — 0, so there exists to > such that 

|7V^^^(t,0)| < ^, ioTt>to. 

Let ( G C[0, 1] and T G C{n) be given. Since the solutions of Equation IT^ are continuous, then there 
exists S > such that, if 

|C(m) -71(0)1 < (5 and |T(m, a) - F(0, a)| < 5, 
for m G [0, 1] and a G [0,Tmax], then 

|iV^'^(t,m)-iV^'i^(t,0)| < ^, 

for t e [to,to + Tmax] and m G [0, 1]. 
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Now, since /i and T are continuous, then there exists b € (0, 1) such that 

|/I(to) -7Z(0)| < (5 and |r(m, a) - r(0, a)| < 5, 
for m e [0,5] and a G [0,T,„a2;]- We define the foUowing functions, 

77 ('.r,^ - / if"^e[0,5], , J r(TO, .), ifme[0,5], 

ifme[&,i], ^^(™")-\ r(6,.), if me [6,1]. 

Then, for m e [0, 1] and a e [0, Tmaa;], we get 

[Mb (m) -71(0)1 < (5 and |rfc(m, a) - r(0, a)| < <5. 

Consequently, 

\NT'^^''^it,m)-NT'-^{t,0)\ < ^, 

for is [ta,to + Tmax] and m G [0, 1]. It foUows that 

\N'^''-^''{t, m)\<v, for te [to, + Tmaa;] and TO e [0, 1]. 

Since N^''-^'' (t, to) is a sohition of Equation (fT^ and since this equation is autonomous for t large enough, 
then N^i"^''{t,m) becomes an initial condition of (fT^ on [to, to + Tmax] x [0,1]- We deduce, from 
that there exists T > to + Tmax such that 

|^'^'"^''(t,TO)| < e, for t > T and m e [0, 1]. 

From Theorem 13.11 there exists t > such that 

N'P'"'^''{t,m) = N'P^^{t,m), for t > t + r^^^ - In (/i(.g(l))) and to G [0,1]. 

Hence, _ 

||7V'^''"(t, .)|| < e, for t > max {T,t + w -In {h{g{l)))}. 



Then, the trivial solution of Equation (jT2l) is globally stable. 

By the same way and using Corollary [312 we show that, if Jl and F are such that 

lim P'^'^(t,0) = 0, 

t — >+oc 

then, the trivial solution of Equation (fT^]) is globally stable. This completes the proof. □ 



Remark 2. One has to notice that the result in Theorem \5.1\ allows us to obtain the global exponential sta- 
bility of the trivial solution of Problem (T^-(W\) on the set Qgs, when this solution is locally exponentially 
stable. The proof in this case, is identical to the previous one. 

The behaviour of the immature cells population has been studied in Section [H The local stability 
of System (fT2|l - (|13p has been studied by Adimy and Crauste [2]. The author proved that, under the 
assumptions that the function x ^ Pim, x) is uniformly bounded and the fmiction x t-^ x/3(to, x) is 
locally Lipschitz continuous for all m G [0, 1], then the trivial solution of p2p - (fT5)) is locally exponentially 
stable if 

1 + 2 sup C{t,m)] sup f /3(to, 0)) < min i inf (s{m) + V'{m)], inf ("f{m) + V'{m)]\. {39) 
(m.t)ef2A / me[o,i] ^ ^ t™e[o.il^ ^ ™e[o,i] V /J 

The proof is based on an induction reasonning. 

Then, we can deduce the following corollary, which deals with the global stability of the system. 
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Corollary 5.1. Assume that Condition ^14^ and Inequality i39\) hold. Then the trivial solution of System 
US ^) - !T^) is globally exponentially stable. 

Proof. From Inequality ([M)) . we obtain the local exponential stability of the trivial solution of (|12p -(fT 
From the definition of p, we get 

inf (s(m) + V'(m)] < p. 

m£[0,l] V / 

Moreover, 



?(0) < sup ^(t,m) and ;9(0,0)< sup /3(m,0). 

(m,t)60A mG[0,l] 

Hence, we obtain 

(2C(0)-l)/3(0,0)<p. 
Then, Theorem 14. II vields to the global stability of the trivial solution of ((^^ - ((50)) 

By using Theorem 15. II and Remark O we conclude. □ 

As an example, let us suppose that d > and 7 > are constant, V and g are linear functions of the 
maturity m, given, for m G [0, 1], by 

V{m) = m, and gim) — — m, with k > 1, 

K 

and the function /3 is a Hill function (see Mackey [HI), defined by 

9"(m) 



P{m,x) = Poim)- 



with pQ and 9 two continuous and positive functions on [0, 1], and n > 1. One can remark that, in this 
case, the function V satisfies Condition ([1]). 

Furthermore, we assume that the function r is given, for m G [0, 1], by 

T(m) — ln(m + a), with a > 1. 
In this case, r is increasing. Therefore, Condition ^ is satisfied. We obtain that 

A(m) = — (\/ Anm + — aV for to G [0, 1], 



2 

and the characteristic curves are given by 

x(s,to) — TOe*, for s < and m G [0, 1]. 

By remarking that A (to) < to for to G (0, 1] if and only if a > k, then we obtain that the trivial solution 
of P^ - (|13p is globally exponentially stable if 

(1 + 2k) sup /3o(to) < min{(5, 7} and a>K. 
me [0,1] 

In the next section, we conclude our asymptotic study by giving a result of instability, based on the 
results of Section [H 
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6 Instability 

The trivial solution of (fT^ is unstable if it is not stable, this means, if there exists e > such that, for all 
v > 0, there exists (/I, F) G C[0, 1] x C{fl) which satisfies 

llAill < I' and ||r||o < I/, 

and 

\\N'P'^it,.)\\> s, fort>0. 

In the next theorem, we show that the instability of the immature cells population leads to the instability 
of the entire population. 

Theorem 6.1. Assume that 

p< (2e(0)-l)/3(0,0). (40) 
Then, the trivial solution of Problem is unstable. 

Proof. From Theorem 14.11 and (00]), we obtain that the trivial solution of Equation ([^H) is unstable. That 
is there exist e > 0, /I G C[0, 1] and F G C{Q) such that N^'^{t, 0) does not tend to zero when t goes to 
infinity. Then, there exist £ > and {tn)neN, with i„ — > +oo, such that 

m'^{t„,0)>e, fornGN. 

Let us suppose, by contradiction, that the trivial solution of ()12p is stable. Then, in particular, there exist 
ly > and T > such that, if 

||7Z|| < ^ and ||F||o < I/, 

then 

\\m'^{t,.)\\<e, fort>T. 

Consequently, 

1^'^'^ (in, 0)1 < e, for n G N such that t„ > T. 

Since we can choose Jl and F as small as necessary, this yields a contradiction. We deduce the instability 
of the trivial solution of (fT2|) . □ 

One can remark that, even if the trivial solution of (fT2|) is unstable, the trivial solution of (|13|) may be 
stable. 
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